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Abstract. We develop several methods, based on the geometric relationship 
between the eigenspaces of a matrix and its adjoint, for determining whether a 
square matrix having distinct eigenvalues is unitarily equivalent to a complex 
symmetric matrix. Equivalently, we characterize those matrices having distinct 
eigenvalues which lie in the unitary orbit of the complex symmetric matrices. 



1. Introduction 

Our aim in this note is to develop simple geometric criteria for determining 
whether a given square matrix T £ M n (C) is unitarily equivalent to a complex 
symmetric matrix (UECSM). To be more specific, a complex symmetric matrix is a 
square matrix T with complex entries such that T = T* (the superscript t denotes 
the transpose operation) and two matrices A, B £ M n (C) are unitarily equivalent 
if there exists a unitary U € M„(C) such that A = U*BU. 

Our primary motivation stems from the emerging theory of complex symmetric 
operators on Hubert space (see [TJ [2l El HI El El HO] , for instance). To be more specific, 
we say that a bounded operator T on a separable complex Hilbert space 7i is a com- 
plex symmetric operator if T = CT*C for some conjugation C (a conjugate-linear, 
isometric involution) on TL. The terminology stems from the fact that the preced- 
ing condition is equivalent to insisting that the operator have a complex symmetric 
matrix representation with respect to some orthonormal basis [3j Sect. 2.4-2.5]. 

From the preceding remarks, we see that the problem of determining whether 
a given matrix is UECSM is equivalent to determining whether that matrix repre- 
sents a complex symmetric operator with respect to some orthonormal basis. From 
another perspective, we may view our main problem as part of a quest to determine 
the structure of the unitary orbit of the set of all complex symmetric matrices. 

Complicating this endeavor, it is well-known that every n x n complex matrix is 
similar to a complex symmetric matrix [9J Thm. 4.4.9] (see also [H Ex. 4] and [3J 
Thm. 2.3]). It follows that similarity invariants, such as the Jordan canonical form, 
are useless when attempting to determine whether a given matrix is UECSM. This 
greatly complicates our work. For instance, one can show that among the matrices 
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all of which belong to the same similarity class, only the fourth is UECSM. In fact, 
prior to the recent advent of Tener's procedure UECSMTest [11] . only a handful of 
matrices were known to be not UECSM. 

In fact, we are partly motivated by Tener's UECSMTest. His procedure is based 
upon the diagonalization of the selfadjoint components A and B in the Cartesian 
decomposition T = A + iB. Although highly effective, it is often difficult to under- 
stand with this method, in simple geometric terms, why a given matrix is UECSM 
or not. In particular, studying the matrices A and B often gives little insight into 
the eigenstructure of T itself. 

In this note, we proceed along a different route. We develop a number of pro- 
cedures, based upon a direct examination of the eigenstructure of T, for testing 
whether T is UECSM or not. To this end, we require that T has distinct eigenval- 
ues - a condition that is satisfied by all matrices outside of a set of Lebesgue measure 
zero in M n (C). On the other hand, Tener's UECSMTest requires that neither A nor 
B have a repeated eigenvalue. In Section[7l we consider several numerical examples 
and establish that neither our test nor UECSMTest subsumes the other. They should 
therefore be viewed as complimentary procedures. 



2. Preliminary Setup 

Let T be a n x n complex matrix having n distinct eigenvalues Ai, A2, . . . , A„ 
and let u\, U2, ■ ■ . , u n denote normalized eigenvectors of T corresponding to the 
eigenvalues Aj. Since 



det(T* - Ail) = det[(T - A*!)*] = det(T - XJ) = 0, 

it follows immediately that T* has the n distinct eigenvalues Ai, A2, . . . , A rl . Let 
Vi, V2, ■ ■ ■ , v n denote normalized eigenvectors of T* corresponding to the eigenvalues 
Ai. Since eigenvectors corresponding to distinct eigenvalues are linearly indepen- 
dent, it follows that both {ui, U2, ... , u n } and {vi, V2, ■ ■ ■ , v n } are bases for C™. 
Based upon the data 



Ul,U2, ■ ■ ■ ,U n \Vi,V2, 



(1) 



we wish to determine if T is unitarily equivalent to a complex symmetric matrix 
(UECSM). Before proceeding, we require a few preliminary lemmas. 

Lemma 1. Under the hypotheses above we have (v,i,i)j) = whenever i ^ j and 
(ui,Vi) ^ 0. 

Proof. If i =/= j, then 

Xi(ui,Vj) = (X l u l ,v j ) = (Tui,Vj) = (ui,T*Vj) = (ui,XjVj) = Xj{u i: Vj) 

whence (tii,Vj) = since Ai 7^ Aj. On the other hand, if (itj,Uj) = for some i, 
then by the preceding (u^, Vj) = for j = 1,2,..., n. Since {v\, t>2, ■ ■ ■ , v n } is a 
basis for C™, it would follow that (ui,x) = for all x 6 C™ whence Ui = 0. This 
contradiction shows that we must have (ui, Vi) ^ for i = 1, 2, . . . , n. □ 

The following lemma allows us to easily express any x S C" in terms of the bases 
{ui,u 2 , . . . ,u n } and {vi,v 2 , ■ ■ ■ ,v n }: 
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Lemma 2. The following formulas hold for all x € C n : 

(2) 



(x,Uj) 



3 



^-^Uj. (3) 



3 ~i (u-j,vj) 

Proof. By symmetry, it suffices to prove ©• Since {u±, «2, . . . ,it n } is a basis for 
C" and since the expression ^ is linear in x, it suffices to verify ([3]) for x = 
u\, U2, . ■ ■ ,u n - Since (ui,Vj) = if i 7^ j and (ui,v,i) 7^ 0, (J3J can be verified 
immediately by setting x = u,. □ 

Lastly, we require a few words about a useful and practical way to view the 
property of being UECSM. 

Definition. A conjugation on C™ is a conjugate-linear operator C : C™ — » C n 
which is both involutive (i.e., C 2 = /) and isometric (i.e., (Cx,Cy) = (y,x) for all 

x, y eC")Q 

In particular, T is a complex symmetric matrix if and only if T is J-symmetric 
(i.e., T = JT*J), where J denotes the canonical conjugation 

J(zi,Z2,...,z n ) = (¥[,z5,...,z^) (4) 

on C™. Moreover, the most general conjugation on C™ is easily seen to be of the 
form C = SJ where S is a complex symmetric unitary matrix. Lastly, it is not 
hard to show that T is UECSM if and only if T is C-symmetric with respect to 
some conjugation C. 

3. The angle test and its relatives 

In this section we briefly outline several convenient geometric conditions which 
are necessary for a given n x n matrix T to be UECSM (unfortunately, none of 
these procedures is sufficient - see Example [S]). Building upon this material, we 
present a condition in Section [5] which is both necessary and sufficient. 

Recall that T is UECSM if and only if there exists a conjugation C on C™ such 
that T = CT*C ' . If this holds, then it follows easily that 

(t - xiy x = o (t* — xiy (Cx) = o. (5) 

Maintaining the notation and conventions of Section [51 we see that if T is C- 
symmetric, then the conjugation C maps the one-dimensional eigenspace of T cor- 
responding to Xi onto the one-dimensional eigenspace of T* corresponding to Aj. 
This is where we invoke the hypothesis that the eigenvalues of T are distinct. Since 
C is isometric and the vectors itj and Vi are normalized, it follows that there are 
unimodular constants aj such that 

Cui = a t Vi 

for i — 1,2, ... ,n. Since C is isometric, this implies that 

(ui,Uj) = (Cuj,Cu t ) 
= (ajVj,a t Vi) 



4n light of the polarization identity, this is equivalent to ||Cai|| = ||x|| for all x S 
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= ajai{vj,Vi) (6) 

for 1 < n. Taking absolute values in the preceding and utilizing symmetry 

yields the following test which can be implemented easily in Mathematical 

Theorem 1 (Angle Test). Suppose that 

(i) T is a n x n matrix with distinct eigenvalues Ai, A2, . . . , A n , 

(ii) 1*1, 1*2, . . • , «n denote normalized eigenvectors of T corresponding to the 
eigenvalues Aj, 

(iii) v±,V2, ■ ■ ■ ,v n denote normalized eigenvectors of T* corresponding to the 
eigenvalues Aj. 

Under these hypotheses, the condition \(ui,Uj}\ = |(«t,Vj)| for all 1 < i < j < n is 
necessary for T to be UECSM. 

In light of the fact that Theorem Q] takes into consideration the (complex) an- 
gles between the eigenspaces of T and compares them to the (complex) angles 
between the eigenspaces of T* , we refer to the procedure introduced in Theorem [T] 
as the Angle Test. One can interpret the Angle Test as asserting that the geometric 
relationship between the eigenspaces of T must precisely mirror the geometric re- 
lationship between the eigenspaces of T* . In some sense, T and T* must be perfect 
mirror images of each other. In Section [5l we present a refined version of Theorem 
[T] which yields a necessary and sufficient condition for T to be UECSM. 

It turns out that the same principles can also be used in certain cases where the 
eigenvalues of T are not distinct. For instance in [H Ex. 7], a similar argument is 
used to show that the matrix 

is not UECSM whenever |o| ^ \b\. 

The condition |6]) can also be interpreted in terms of Gram matrices. Let 
U = (ui\u2\ ■ ■ ■ \u n ) and V = (i>i|i>2| • • • \v n ) and observe that ([6]) is equivalent 
to asserting that 

{U*Uf = A*(V*V)A (7) 
holds where A — diag(ai, 012, ... , a n ) denotes the diagonal unitary matrix having 
the unimodular constants a\, 0.2, ■ ■ ■ , otn along the main diagonal. This leads us to 
the following test: 

Corollary 1 (Grammian Test). A necessary condition for T to be UECSM is that 
U*U and V*V have the same eigenvalues, repeated according to multiplicity. 

Proof. If T is UECSM, then © holds. Since U*U is a positive matrix, it follows 
that U*U and {U*U) t are both unitarily equivalent to the same diagonal matrix, 
whence ([7]) implies that U*U and V*V are unitarily equivalent. □ 

We should remark that Example [5] in Section [6] reveals that passing the Gram- 
mian Test is insufficient for a matrix to be UECSM. On the other hand, we show 
in Section [4] that ([7|) is both necessary and sufficient for T to be UECSM. 

Taking the determinant of both sides of ([7]) immediately yields the following: 

Corollary 2 (Parallelepiped Test). Maintaining the notation above, if \detU\ 7^ 
|detV|, thenT is not UECSM. 
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The name Parallelepiped Test stems from the fact that | det C/ 1 i and |dety|i 
can be interpreted as the volumes of the generalized parallelepipeds in C™ spanned 
by the vectors u\, U2, ■ ■ ■ -,u n and v%, ■ ■ ■ , v n , respectively. 

The following example illustrates the preceding ideas: 

Example 1. We claim that the matrix 



is not UECSM. Letting Ai 
normalized eigenvectors 



Ml = 



T 



0, A 2 




1, and Xr 



U2 




and 



of T and T* , respectively. Setting 



U 



we immediately find that 





1 




(: 












1 

VE 



2 



V = 










1 






Vl = 










3 















_ 2 
2" 

I 

3 



(9) 



det 17 1 



det V I 



whence it follows from the Parallelepiped Test that T is not UECSM. 
Moreover, we also have 



(u*uy = 



( 1 


1 


1 




? 


1 


1 










1 


1 



v*v 




(10) 



whence, by considering the moduli of the off-diagonal entries in (fTD) , it is clear that 
no diagonal unitary matrix A exists which satisfies (J7J. Thus the Grammian Test 
once again establishes that T is not UECSM. 

Finally, let us take this opportunity to illustrate the Angle Test, which is less 
computationally intensive than either the Parallelepiped Test or the Grammian 
Test. A short calculation based upon the data © and (0 reveals that 

1 2 

\(ui,u 2 )\ = -j= ^ - = |(Ul,f2)| 

whence T is not UECSM. This can be also seen directly by examining the (1,2) 
entry of the matrices in (HI 
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It is important to remark that none of the conditions described in Theorem [TJ 
Corollary[l] or Corollary[2j are sufficient for T to be UECSM. This is illustrated in a 
series of rather involved computations (see Example[5]) that we postpone until later. 
In Section [5] we remedy this situation and provide a test which is both necessary 
and sufficient. 

4. Constructing a Conjugation 

Under our running hypotheses, it turns out that the condition ([6]) is sufficient for 
T to be UECSM (in particular, so is the Gram matrix condition ([7])). The following 
lemma is the main workhorse upon which the rest of this note is based: 

Lemma 3. Let 

(i) T be a n x n matrix with distinct eigenvalues Ai, A2, . . . , X n , 

(ii) 1*2, . . • , «n denote normalized eigenvectors of T corresponding to the 
eigenvalues Xi, 

(iii) v±,V2, ■ ■ ■ ,v n denote normalized eigenvectors of T* corresponding to the 
eigenvalues Xi. 

If unimodular constants ai, 0:2, . . . , a n exist such that 

(ui,Uj) =a~iCij{vj,Vi) (11) 

holds for\<i <j < n, then T is UECSM. 

Proof. First observe that if (fTTj) holds for 1 < i < j < n, then (JTTJ) holds whenever 
1 < i,j < n by symmetry and the fact that we are considering normalized eigen- 
vectors. Let Cui = onvi for i — 1, 2, . . . , n and extend this by conjugate-linearity 
to all of C™. We intend to show that C is a conjugation with respect to which T is 
C-symmetric. Since C is conjugate-linear by definition, it suffices to show that C 
is involutive and isometric. 

Step 1: Show that C is involutive (i.e., C 2 = I). 

Since C 2 is linear, it suffices to verify that C 2 Ui — m for i = 1, 2, . . . , n. By © 
it follows that 



n 1 1 

Cui = a l v l = ai - - Uj (12) 



whence 



c 2 Ul = c \ ai y pi^J-uA bym 



CCi} 7—^ — % \cuj conjugate-linearity 

E(Vj,Vi) 
rCt-iVj definition of C 

Eaiouiui.Uq) , , 

h ,\ ) a > Vj by m 
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= Ui. by |HJ 

Thus C is involutivc. 

Step 2: Show that C is isometric (i.e., ||Ca;|| = ||a;|| for all x G C"). 
If x = ^27=1 c i u ii then observe that 
\\x\\ 3 = {x,x) 

n 
n 

= ^2 CiCjcHaj(vj,Vi) by dXTJ) 

i,J=l 
n 



CiCj{ajVj,a.iVi) 



= CiCj(Cuj,Cui) definition of C 



(cjCuj,CiCui) 



i=l i=l 



= (Cx, Css) 

= l|c^ll 2 - 

Thus C is isometric whence C is a conjugation on C™. 

Step 3: Show that T is C-symmetric (i.e., T = CT*C). 

Since both T and CT*C are linear, it suffices to prove that they agree on the 
basis Ui, U2, ■ ■ . , u n . Having shown that C 2 = I, it now follows from the equation 
Cui = ctiVi and the conjugate-linearity of C that Cvi = OfjUj. Thus 

CT*Cui = CT* '{ctiVi) — WiCT*Vi — alCXiVi 
= aiXiCvi = aiXiCtiUi = A* it* = Tm 

whence T is C-symmetric and hence UECSM. □ 

The conjugation C constructed by Lemma |3] can be concretely realized as C = 
SJ where S is a complex symmetric unitary matrix and J denotes the canonical 
conjugation ([4} on C". Let us briefly describe the construction of the matrix S. 

First observe that C satisfies Cui = c^Vi, which is easily seen to be equiva- 
lent to Cvi = ctiUi for i = 1,2, ...,n. As before, let U = (ui |t*2 1 ■ ■ ■ \u n ) an d V = 
(ui|u2| • • • \v n ) denote the matrices having the vectors u\, ■ ■ ■ , u n and vi, «2, • • • , v n , 
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respectively, as columns. Since the columns of U and V form bases of C™, it follows 
that both of these matrices are invertiblc. Next we note that 



v*u 



( (wi,«i> 



V 2 



{ui\u 2 \ ■ ■ ■ \u n ) 



by Lemma [1] Let 



{U 2 ,V 2 ) 



\ 



E 



(ot x 



Q'2 



(13) 



and 



D = AE' 1 



1 



V 



Q2 
(«2,U2) 



(«".«») / 



We claim that C = S 1 J where 



5 = UDU 



(14) 



(15) 



To prove (fl~5|) , it suffices to show that the conjugate-linear operators C and S J 
agree on each of the vectors Vi. In other words, we must show that SJvi = aiui 
for i = 1,2, ... ,n. Letting si, S2, . . . , s n denote the standard basis for C™ we have 

SJvi = UDU l Jvi 

= UAE-Wjvi 

= VAE- x J\J*v, 



= UAJE 

= UAJV^ 1 
= UAJs r 
= UA Sl 
= UotiSi 

= CXiUi. 



U* Vi 



by HI]) 

by (O 

since JU* = U t J 
l 



since JE 

since V^ 1 = E 
def. of V 
since Jsi = Sj 
by ([13]) 
def. of [/ 



E- X J 



U* 



Thus C = SJ . Since the matrix D = AE -1 is diagonal, it is clear from (p~5|) that 5 
is symmetric. Since S* = CJ is the product of two conjugations, it is an invertible 
isometry and hence unitary (see also Lem. 1]). 

It is worth remarking that the condition T = CT*C implies that T = SJT*SJ = 
STtJSJ = ST'S 1 * since S is symmetric (i.e., S is J-symmetric). Therefore the 
matrix S yields a unitary equivalence between T and its transpose T*. 



5. The Strong Angle Test 



The main theorem of this article is the following necessary and sufficient condi- 
tion for a matrix with distinct eigenvalues to be UECSM. The procedure introduced 
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in the following theorem was implemented in Mathematica by the first author. We 
refer to this procedure as StrongAngleTest. 

Theorem 2 (Strong Angle Test). If 

(i) T is a n x n matrix with distinct eigenvalues Ai, A2, . . . , A n , 

(ii) tii, u>2, ■ ■ • , u n denote normalized eigenvectors of T corresponding to the 
eigenvalues Ai, A2, . . . , X n , 

(iii) Vi , v%, . . . , v n denote normalized eigenvectors of T* corresponding to the 
eigenvalues Ai, A2, . . . , X n , 

then T is UECSM if and only if the condition 

(Ui,Uj)(Uj,U k )(Uk,Ui 



{v il v j )(v J1 v k ){vk 1 v l ) 
holds whenever l<i<j<k<n and not all ofi,j,k are equal^ 



(16) 



Proof. The necessity of the condition (|16p follows immediately from ([6]) . The proof 
that (fTr?)) is sufficient for T to be UECSM is more complicated. First observe that 
if (fl6|) holds for l<i<j<k<n, then (fTB"|) holds whenever 1 < k < n by 



symmetry. Let us assume for the moment that (u^Uj) ^ (whence (Vj,Vi) 7^ 0) 
for 1 < i,j < n. Later we will relax this restriction, but for the sake of clarity it 
will be easier to consider this special case first. Under this additional hypothesis, 
there exist n 2 unimodular constants 0ij uniquely determined by 



{Uj,Uj) 



(17) 



for 1 < i,j < n. Since = \\vi\\ — 1 by hypotheses (ii) and (iii), it follows 
immediately that that /3a = 1 for 1 < i < n. Moreover, we also have 



0ij{Vj,Vi) = (Ui,Uj) = (Uj,Ui) =0ji{Vi,Vj) = [3ji{Vj,Vi), 

whence /3y = In other words, the matrix B = (/3y )™ J=1 is selfadjoint and has 
constant diagonal 1. Suppose for the moment that, based on the hypothesis (|16[) . 
we are able to establish that B enjoys a factorization of the form 

/ 1 012 013 ■ : X 

021 1 023 ■ 
031 032 1 ' 



02n 
03n 



\0nl 0n2 0n3 





( al \ 




02 




\a n J 



OCn) 



(18) 



(i.e., suppose that we are able to show that B is positive and has rank one). By 
(jTT)) and the preceding factorization (TTHJ) it would then follow that the unimodular 
constants 011, a.2, ■ ■ ■ > c^n satisfy 

(ui,Uj) =alaj(vj,Vi) (19) 

for 1 < i,j < n. At this point, we could invoke Lemma [3] to conclude that T is 
UECSM. 



2 Observe that setting k = j in condition I I16H leads to \(ui,Uj)\ = \(vi,Vj)\ for 1 < i < j < n. 
Thus Theorem [2] can be viewed as an extension of the original Angle Test (Theorem [TJ . Also 



note that if i = j = k, then I I16H merely asserts that \\ui\\ 
conditions (ii) and (iii). 



which is already known from 
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The difficulty in the approach outlined above lies in the fact that some of the 
inner products (v,i,Uj) or (vj,Vi) may vanish. If this occurs, then we cannot im- 
mediately consider the associated unimodular constants 0ij defined by (jTTJ) since 
applying (|16|) with k = i implies that (vj, Uj) = if and only if (itj, Uj) — 0. On the 
other hand, observe that the hypothesis (fl6|) implies that 



a a a (uj,U k ) (u k ,Uj) 

{Vj,Vi) {Vk,Vj) {Vi,v k ) 

holds whenever 0ij, 0j k , k i are well-defined by (jTTJ) . In light of the fact that each 
flij is unimodular, we obtain the following multiplicative property 

/3y = PikPkj (20) 

whenever the expressions above are well-defined by (|17p . 

Regarding the matrix B = (0ij)fj—i as being only partially defined by (JTTJ), 
suppose for the moment that we are able to define unimodular constants 0ij for 
those i and j for which (itj, u,-) = (u,-, Uj) = such that the multiplicative property 
(f2T))) holds for all 1 < i, j, k < n. Under this hypothesis, we claim that the matrix 
B = (flij)?j=i has a factorization of the form (fT5|) . Indeed, use (|20p and the fact 
that (3ji = Pij for 1 < i, j < n to conclude that 



B = 



( 1 

021 



012 
1 

/3 32 



013 
023 
1 



03n 



0V. 



0: 



13 



(1 012 /? 



13 



/?!„) ■ (21) 



\0nl 0n2 Pn3 ■■■ 1 / \01 n J 

As suggested by (|18p . we now define the unimodular constants ai,a2, 
setting ai = 0u for 1 < i < n. Next observe that 



,a n by 



aittj = 0u0\j = 0n0ij = 0i_ 



(Ui,Uj) 



Uj, 



holds whenever 0ij is well-defined by (fl"7|) . Thus the desired condition (|19[) holds for 
all 1 < i,j < w (since it holds trivially if (Ui,Uj) — (vj,Vi) — 0) and T is UECSM 
by Lemma [3] 

To complete the proof of Theorem [2] it suffices to demonstrate a procedure by 
which we may define unimodular constants 0ij for those i and j for which (m, Uj) — 
(vj,Vi) = such that the multiplicative property ([2H)l holds for all 1 < i,j < n. 



This will lead us to the desired matrix factorization (|2Tj) . 

To define the constants 0^ we employ an inductive procedure. Consider the 
partially defined n x n matrix 



1 


012 


013 ■ 


• 01r 


* 


* 


021 


1 


023 ■ 


■ 02r 


* 


* 


031 


032 


1 


• 03r 


* 


* 


Prl 


0r2 


0r3 ■ 


1 


* 


* 


* 


* 


* 


* 


1 ■ 


* 


* 


* 


* 


* 


* 


■ 1 



\ 



(22) 
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where * indicates either an entry (3ij already defined by (|17j) or an entry that is not 
defined in terms of (fT7|) because (ui, Uj) — (vj,Vi) — 0. As our inductive hypothesis, 
we assume that the multiplicative property (|20p is satisfied by all triples /3y , 0ik , 0kj 
for which 1 < i, j, k < r. 

To complete the proof Theorem [2] we must devise a way to fill out the undefined 
entries in ([22| with unimodular constants in such a way that (|20|) holds for 
these new entries. There are two cases to consider: 

Case 1: Suppose that there exists an entry /3 i( v +1 ) with 1 < i < r in (|22]) that 
is already defined by (fT?)) . Without loss of generality, we may assume that it 
is the /3i( r -fi) is the entry that is well-defined by (|17p since this situation may 
be obtained by permuting the indices 1,2, ...,r and relabeling the eigenvectors 
U\, U2, ■ ■ ■ ,u r ; v\, t>2, ■ ■ ■ , v r . We are thus left with the partially completed matrix 

/ 1 012 Pl3 ■■■ Plr /3l(r+l) * ' • ' * \ 



1 


/3l2 


fa ■ 


• Plr 


Pl(r+1) 


* 


* 




1 


P23 ■ 


■ 02r 


* 


* 


* 


031 


032 


1 


■ 03r 


* 


* 


* 


Prl 


A-2 


Pr3 ■ 


1 


* 


* 


* 


/3(r+l)l 


* 


* 


* 


1 


* 


* 


* 


* 


* 


* 


* 


1 ■ 


* 


* 


* 


* 


* 


* 


* 


■ 1 / 



V 

For each entry A( r +i) with 2 < i < r (i.e., the entries immediately below /3i( r +i) 
and above the 1 on the main diagonal) there are two possibilities: 

Subcase 1.1: If /3i( r +i) is already well-defined by (fl7)) . then do nothing. 

Subcase 1.2: If Pit r +i) cannot be defined by (fl"7|) because (ui,Uj) = (vj,Vi) = 0, 
then let 



Pi(r+i) PaPi(r+i) 
to obtain the partially defined matris0 

/ 



(23) 



1 


012 


013 


Plr 


Pl(r-{ 


-1) 


* 


* 


#21 


1 


023 


02r 


02(rA 




* 


* 


031 


032 


1 


03r 


03(rA 


'1) 


* 


* 


Prl 


0r2 


0r3 


1 


0r(r-\ 


-1) 


* 


* 


P(r+l)l 


/3(r+l)2 


0{r + l)3 


0(r + l)r 


1 


* 


* 


* 


* 


* 


* 


* 




1 ■ 


* 


* 




* 


* 


* 




* 


• 1 ( 



(24) 



Case 2: Suppose that there does not exist an entry 0i( r +i) with 1 < i < r in 
([22]) that is already defined by (fT7|) . In other words, suppose that (m 7 u r+ i) = 



^The entries /3( r _|-i)i, P( r +l)2> /3(r+l)3i • • • : Pr(r+V) in the (r + l)st row are denned by conjugate 
symmetry: ftj = (3ji- 
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(v r +l,Vi) 

matrix 



= whenever 1 < i < r. We are thus left with the partially completed 

/ 1 012 013 • ■ ■ 01r * * • ■ ■ * \ 



V 



1 


012 


013 ■ 


■ 01r 


* 


* 


* 


021 


1 


023 ■ 


■ 02r 


* 


* 


* 


031 


032 


1 


■ 03r 


* 


* 


* 


0rl 


0r2 


0r3 ■ 


1 


* 


* 


* 


* 


* 


* 


* 


1 


* 


* 


* 


* 


* 


* 


* 


1 ■ 


* 


* 




* 


* 


* 


* 


■ 1 



/ 



Select a complex number of unit modulus and assign this value to Having 
done this, we define @i( r -\-\) for 1 < i < r as in (|23[) to obtain a partially defined 
matrix of the form (1241). 



To wrap-up the proof, we must show that in either case (|23[) defines the new 
entries /3 i ( r+1 ) in a manner which is consistent with the multiplicative property 
(f2"TJ)) . For 1 < i, k < r we employ the definition (f2"3"|) to find that 



/3ikPk(r+i) = {Pn0ik)(0kiPi(r+i)) b y inductive hypothesis and (|23J) 

= Pa ifi\kPxk)P\(r+i) hermitian symmetry 

= PnPi( r +i) since \(3i k \ = 1 

= A(r+i)- by ([231) 

Thus (f23|) defines /3;( r+ i) for 1 < i, j < r + 1 in a manner consistent with (|20|) . 

Starting with the upper left lxl block, repeated applications of the preceding 
inductive procedure eventually yields an n x n matrix B = (/%)™ 7 - = i whose en- 



tries satisfy the required multiplicative condition 
Theorem [2j 



This concludes the proof of 

□ 



Corollary 3. Every 2x2 matrix is UECSM. 

Proof. Let T be a 2 x 2 matrix. If T has a repeated eigenvalue A, then by Schur's 
Theorem on Unitary Upper Triangularization, it follows that T — XI is unitarily 
equivalent to a scalar multiple of a 2 x 2 nilpotent Jordan matrix. This Jordan 
matrix is C-symmetric with respect to C{z\,Z2) — (~Z2,~zi) whence T is UECSM. 
We therefore restrict our attention to the case where T has two distinct eigenvalues. 
Upon applying Schur's Theorem, subtracting a suitable multiple of the identity, and 
normalizing, we may assume that 

H <25) 

for some complex constant a. A short computation reveals that normalized eigen- 
vectors for T corresponding to the eigenvalues Ai = and A2 = 1 are 
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Similarly, we find that corresponding normalized eigenvectors for T* are given by 

t* = NT ) , * = (?) ■ (27) 

By Theorem El T is UECSM if and only if fijfy holds for all 1 < i < j < k < 2 
such that not all of i,j, k are equal. This leaves us only two ordered triples k) 
to consider: (1, 1, 2) and (1, 2, 2). These values of k both lead to the condition 
U2)| = 1(^1,^2)1 which needs to be verified. Since 

1(^2)1= H =l(ttl,m)l 

v 1 + M 

follows immediately from ((26]) and (EZj), we conclude that T is UECSM. □ 

The preceding corollary has been proved in a number of different ways by several 
different authors. For instance, one can reduce to the special case (|23|) as above and 
then construct the corresponding conjugation by straightforward computation [H 
Ex. 6]. The procedure developed by J. Tener can also be used to establish Corol- 
lary [3] [HJ Cor. 3]. We should also mention that Corollary [3] is the byproduct of 
more sophisticated theorems. For instance, it follows immediately from N. Chevrot, 
E. Fricain, and D. Timotin's study of the characteristic functions of complex sym- 
metric contractions [TJ Cor. 3.3]. More recently, the second author and W. Wogen 
established that every binormal operator (i.e., an operator that is unitarily equiv- 
alent to a 2 x 2 block operator whose entries are commuting normal operators) is 
complex symmetric [6]. Corollary [3] is a special case of this result. 

6. A FEW EXAMPLES 

To illustrate the preceding ideas, we devote this section to the detailed consider- 
ation of several examples. In particular, Example [5] demonstrates that none of the 
simple conditions given in Section [3] is sufficient for T to be UECSM. 

Example 2. Let T be a 2 x 2 matrix with eigenvalues Ai 7^ A2. As before, let 
Hi, U2 denote normalized eigenvectors of T corresponding to the eigenvalues Ai, A2, 
respectively and let v\ , V2 denote normalized eigenvectors of T* corresponding to 
the eigenvalues Ai,A2, respectively. 

By Corollary^ we know that T is UECSM and hence \ (ui,u 2 } \ = |(u 2 ,i>i)|. We 
may therefore find unimodular constants ot\ and c*2 such that 

(ui,u 2 ) = aTa2(v2,vi). 

For instance, if T is normal, then we may simply set ot\ = 02 = 1 since (111,112) — 
(v2,vi) — 0. Letting 



U = 



Mil U12 
U21 U 2 2 

denote the matrix whose columns are the eigenvectors 1*1,1*2 of T we follow the 
procedure outlined at the end of Section 3] to construct the conjugation 



C 



I OL\u\ x 02^12 O^UuVai ^ 02_Ul2^22 \ 

(Ul,l>l) (U2,V 2 ) (ui,Vl) (U 2 ,V 2 

CK1U11U21 a 2 Ml2«22 Cellini ! «2W| 2 



(28) 



\ (tll,Ul) {u 2 ,v 2 ) (tll,Vl) (u 2 ,v 2 ) I 
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with respect to which T is C-symmetric. 

Example 3. Applying the preceding formula to the matrix (|25p . where a ^ 0, and 
using the data ([26]) and (|27|) we find that 



{111,112) 



Following the notation of Example note that 

"11 = 1 



u 2 i = 



U12 



u 2 2 



1 



One possible solutiorQ to the equation 

(ui,u 2 ) ^cna 2 (v2,v 1 ) 

is given by 

ai = FT) "2 



a a 
\a\ \a\' 
Plugging this data into P5)l we obtain the conjugation 

H_ 

' x ' 



C 



/ a/\a\ 
[oj 



-g/N 



with respect to which the matrix (|25|) is C-symmetric. 
Example 4. In [Til Ex. 3], the matrix 

/0 7 
T = I 1 -5 
\0 6 

is demonstrated to be UECSM via Tener's UECSMTest. For the sake of comparison, 
let us also consider this matrix using the techniques discussed above. Letting Ao = 
6, Ai = 1, and A2 = 0, we obtain the corresponding normalized eigenvectors 



"1 



of T and 



vi 




U3 = 



f'3 




of T* , respectively. A short computation reveals that 

(V1,V 2 ) = -7=, (V2,V 3 ) 



(Ul,U 2 ) 



5V2' 
7 



y/2' 



5V? 



(V3,vi) 



7 

IT 



7 

IT' 



^The other solutions will simply yield a unimodular multiple of S. 
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whence it is clear that (|16|) holds for all triples 

k) = (1, 1, 2), (1, 1, 3), (1, 2, 2), (1, 2, 3), (1, 3, 3), (2, 2, 3), (2, 3, 3) 

required by StrongAngleTest. In particular, this proves that T is UECSM. 
The corresponding matrices U — (ui|u2| u 3) an d V = (vi^lv^) are 



U = 





7 














\ 11 





















1 








V2 



42 
55 
7_ 
11 



As expected, T passes the Parallelepiped Test (Corollary [2]) since 

3V2 



Idetm = I det VI = 



55 



Next, observe that 



U*U = 



l 



l 

V2 



l —4b --h 



i -4 4t 



l 



V*v= -fe 



J7_ 

11 



7 

5V2 



1 



J7_ 

11 





ai 






















Q!3 



whence T passes the Grammian Test (Corollary [TJ with the A from |(7J) being 



In particular, this once again confirms that T is UECSM. 

Let us, for the moment, examine the mechanics of the proof of Theorem^ which 
establishes the theoretical underpinnings of the procedure StrongAngleTest. Using 
the data above, we find that the matrix B = (/3y ) from the proof of Theorem [2] is 
given by 



B 



(u*uy 
v*v 



1 



1 




1) 



whence we again read that a\ — l,ct2 — — l,aa = — 1. We remind the reader 
that the quotient appearing in the preceding equation is simply the entry-by-entry 
quotient of the matrices (U*UY and V*V. 

Based upon the preceding calculations, we can construct the corresponding con- 
jugation C = SJ where S is a complex symmetric unitary matrix which is to be 
determined (this matrix also has the property that T = ST f S*). Following the 
recipe described at the end of Section 2] we obtain 



E = V*U = 



so that 



D = AE~ 



6 
11 










1 

10 


i) 





- 


55/ 


6 





0\ 





-10 




V« 
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Putting this all together we find that 




S = UDU* = 



It turns out that our S differs from the corresponding matrix obtained in Ex. 3] 
by a unimodular multiplicative factor of ~ i9+6z\/74 ^ 

The following important example demonstrates that the Angle Test (Theorem 
[lj, the Grammian Test (Corollary QJ, and the Parallelepiped Test (Corollary [2|) 
are insufficient to determine whether a given matrix is UECSM. In particular, this 
demonstrates the utility of the Strong Angle Test (Theorem [2|) , which provides a 
necessary and sufficient condition. 

Example 5. Consider the matrix 



/5 


-1 


3 \ 


2 4 


f 


2 


2 -2 


6 


-2 


\0 -2 


f 


4/ 



which has the distinct eigenvalues 

Ai = 5 + iV§, A2 = 5 — i\/5, 
and corresponding unit eigenvectors 



A 3 = ±(9 + n/l5), A 4 = |(9-n/l5), 



/ 3 3N/5 » 

4 f— 5* - 2^5) 



15 



V 



7! 



/ 2i 1 

' 3 37S 

^(5t-2>/5) 
£(-5i->/5) 
75 



/ »(-5t+Vl5) \ 

2733 


75 ~^ 73o 



V 



/ i(5i+\/l5) \ 
2V30 



75 73o 

v v/S 



The matrix T* has the eigenvalues 

A7=5-iV5, A^=5 + i\/5, A3" = |(9 - n/l5), Ai = |(9 + «Vl5) 
and corresponding unit eigenvectors 




10 



[5i 



2 V / 15 



^(-lSz + v 7 ^) 



V 



/15 



The matrices U = (tti|u2|us|u4) and V = (vi \v2\v3\v4) satisfy 

|detC/| = ^ = |detF| 
whence T passes the Parallelepiped Test (Corollary [5]). 
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A further computation reveals that the matrices 
/ 



U*U = 



4+y/3 
-4+n/3 

6V2 



v*v = 



3V5 

(l + 4v/3") 

6VT0 
(-l+4>/3) 



" 3^ 

i(3 + 4\/3") 



3 3^5 
1 

-4 + V3 _ '(-1 + 4^) 
6V^ 6VT0 

4+V3 _ '( 1 + 4 ^ ? ) 
6\/2 6VT0 

1 _ gj 

3 3^5 

1 



4 + V3 _ 

6V2 6vTo 

(-1 + 4V3) 
6VT0 



-4 + n/3 

6\/2 



-4+y3 _ i(-l+4V3) \ 
6V2 6\/To 

6v% 



6\/T0 



2\/3 6\/5 

i(-3 + 4x/3") 



2V5 



6V5 



2^3 
1 

2V3 ' 



i(-3+4v/5) 

i(3+4\/3) 
6\/5 



2^ 
1 

" 2s/3 



i(3+4v^) 
i(-3 + 4 % /3") 



1 

" 2^3 
1 

2s/3 ' 



;(-3+4v / 3") \ 

i(3+4v/3) 
675 



share the eigenvalues (given approximately by) 

2.73115, 0.932497, 0.253856, 0.0824931. 

Thus T passes the Grammian Test (Corollary [TJ. 

Recall that the ijth entries of U*U and V*V are (uj, Uj) and (vj,Vi), respectively. 
Therefore to check whether T passes the Angle Test (Theorem [1} , we need only 
compare the moduli of the entries of U*U and V*V. The moduli of the entries of 
U*U and V*V are equal, entry- by-entry, and given by 



/ 



^ ( 3 + v^) ^ (3 - 



3+V3 
3 - 

1 

2 
r> 



3-v^ 

3 + ^' 

2 
5 



Thus T passes the Angle Test (Theorem [1}. 
On the other hand, since 

(ui,U2)(W2.«3>(W3)«l) = ^(S-iV 7 ^) ^ ^(5 + iV5) = («l,U2)(«2,f3>(V3,«l), 

the Strong Angle Test (Theorem^ asserts that T is not UECSM. Similar compu- 
tations reveal that the desired condition (fTBf is violated for the triples (i,j,k) = 
(1,2, 4), (1,3, 4), (2, 3, 4) as well. 

Working through the mechanics of the proof of Theorem [2j we find that the 
matrix B = (/3y ) is given by 

/ 



B = 



ll + 6v / 3 + 2iv/5 + 3iv / T5 

4 % /2"(3 + v^i) 
ll-6V3 + 2iv^5-3iy / r5 
\ 4v^(-3+%/3") 



ll-6%/3-2iy'SS + 3ivT5 

4v^(-3 + V3) 
1 1 + 6 \/3- 2iV*5-3iyT5 
4V5(3+x/3) 



ll + 6v^3-2iv / 5-3i^L5 

4v'2(3+v / 3) 
ll-6^/3"+2^v / 5-3z^/l5 
4v / 2(-3 + v / 3") 
1 



ll-6V'3"-2iv / 5 + 3t\/T5 \ 

4n/2(-3+v^) 
ll + 6V3 + 2iV5 + 3ivT5 
4V5(3+%/3") 



7 + IV15 ) 



In particular, each entry of B is unimodular whence we once again see that T 
passes the Angle Test. Also observe that the rank of B is 4 and its eigenvalues are 
approximately 

3.88114, 0.694237, -0.66798, 0.0926015. 

In particular, B is neither rank-one nor positive. 

We should also mention that J. Tener's procedure UECSMTest also confirms, via 
entirely different methods (see Section [7]), that T is not UECSM. 
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The preceding example was discovered by the first author during a search of 10 
million random integer matrices. Such examples appear to be exceedingly rare and 
those which can be worked through in closed form rarer still. Moreover, we were 
unable to find a 3 x 3 matrix with the same properties. 

7. Comparison with Tener's UECSMTest 

J. Tener's procedure UECSMTest, introduced in [TT], is an effective tool in deter- 
mining whether a given matrix is UECSM. However, there are certain limitations 
inherent in the procedure. To be more specific, UECSMTest cannot be applied if the 
given matrix T is 4 x 4 or larger and either Cartesian component A or B in the 
decomposition T = A + iB (where A = A* and B = B*) has a repeated eigenvalue. 

On the other hand, the criterion for applying StrongAngleTest is simply that the 
matrix T have distinct eigenvalues. In this section, we compare the two procedures 
and demonstrate the existence of matrices, both UECSM and not, for which either 
UECSMTest or StrongAngleTest (possibly both) fail to apply. In particular, this 
demonstrates that StrongAngleTest and UECSMTest are complimentary procedures 
in the sense that neither test subsumes the other. 

Obviously many normal matrices (e.g., the 4x4 identity matrix) do not satisfy 
the hypotheses of cither test. This does not pose a problem, however, since the 
Spectral Theorem asserts that every normal matrix is unitarily equivalent to a 
diagonal matrix whence every normal matrix is UECSM. In light of the preceding 
remarks, we therefore focus our attention on producing examples which are non- 
normal. 

Finding non-normal matrices for which StrongAngleTest is applicable and for 
which UECSMTest is not is relatively straightforward. Several examples are listed in 
Table 1 below (where c(T), cr{A), &(B) denote the spectra of the operators T, A, B, 
respectively, in the decomposition T = A + iB, A = A* , B = B*). 



T 


a(T) 


a(A) 


a(B) 


UECSM? 




'2 N 
4 
8 4 

v -2, 


\ 
/ 


-2,2,4,8 


2,4,3± v/29 


0,0, ±2 


Yes 




/2 0\ 
4 
2 

\0 8 0/ 




2,4, -2± 2iV3 


distinct 


0,0, ±V21 


No 




/ 4 1-1 -2\ 
3 2-41 
-1-2 4 1 
V-4 1 3 2 J 




-2,2,4,8 


2,4,3± V29 


0,0, ±2 


Yes 




/ 4 -1 1 -2\ 
-2 1-14 
-14-2 1 

\ 1 -2 4 -l) 




2,4, -2±2«V3 


distinct 


0,0, ±V21 


No 



Table 1: Examples of simple matrices for which StrongAngleTest is ap- 
plicable and UECSMTest is not. The third and fourth matrices listed are, 
respectively, unitarily equivalent to the first and second matrices. The eigen- 
values of the second and fourth matrices are distinct, but too long to display 
explicitly in the confines of the tabic. 
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In cases where T has repeated eigenvalues, one frequently finds that both A and 
B both have distinct eigenvalues. Such matrices are testable by UECSMTest but not 
by StrongAngleTest. Several simple examples are listed in Table 2 below. 



T 




a(T) 



0,0,0 



0,±9V2 



UECSM? 



Yes 



18 





0,0,0 



0,± 



»/5 



0,± 



9v/5 



No 



' 8 + 


4i 


H 


-8i 


-8 


+ 8i\ 


-8- 


{- 2i 


-4 


+ 4t 


8H 


-4j 


v 4- 


Ai 


2 - 


- 8t 


-4 


- Si) 




2i 


4 + 


4i 


-8- 


I-4A 


(H 


+ i 


-4- 


h 2i 


8 + 


2i 




2i 


2- 


4i 


-4- 


- 4i/ 



0,0,0 



0,^:9^2 



0,±9\/2 



Yes 



0,0,0 



0,± 



0,± 



9^5 



No 



Table 2: Matrices for which UECSMTest is applicable and StrongAngleTest 
is not. The third and fourth matrices listed are, respectively, unitarily equiv- 
alent to the first and second matrices. 



It is possible to construct matrices for which neither StrongAngleTest nor 
Tener's UECSMTest is applicable. To be more specific, we exhibit several matri- 
ces T such that 

(i) T has repeated eigenvalues (so that StrongAngleTest is not applicable), 

(ii) T = A + iB is at least 4x4 and either A or B has repeated eigevalues (so 
that UECSMTest is not applicable). 

Although it is relatively straightforward to produce matrices T satisfying (i) and 
(ii), it is naturally quite difficult to check whether T is UECSM or not since by 
design neither StrongAngleTest nor UECSMTest are applicable. Fortunately, the 
set of matrices having properties (i) and (ii) has Lebesgue measure zero in M n (C) . 

We require a couple preliminary lemmas. The following can be found in [111 
Ex. 1] or [1 Ex. 1]: 

Lemma 4. The matrix 

a 0\ 

6 (29) 
0/ 

is UECSM if and only if ab = or \a\ = \b\. 

In particular, the matrix (|29[) is not UECSM whenever a and b are nonzero and 
satisfy \a\ ^ In our construction, we intend to use ([2^)) as a building block in 
conjunction with the following lemma from [7j: 

Lemma 5. T is UECSM if and only i/TffiO is UECSM. 

In the lemma above, T© denotes the orthogonal direct sum of T with a square 
zero matrix of any given size. Since T is UECSM if and only if T — XI is UECSM 
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for any A € C, it follows from Lemma |4] and Lemma [5] that the matrix 



T 






V o 



\ 





b 
o / 



(30) 



can be made UECSM or not according to our choice of a and b (the value of c is 
irrelevant). A short computation then reveals that cr(T) = {0, 0, 0, c} and 



a(A) 
a(B) 



{0,Rec,±VM 2 + |&| 2 }, 
{0,Imc,±v/|a| 2 + |6| 2 }, 



whence if c is either real or purely imaginary condition (ii) holds. This leads us to 
the examples listed in Table 3 below: 



T 


a(T) 


a(A) 


a(B) 


UECSM? 




/4 0\ 
8 
8 

Vo 0/ 




0,0,0,4 




0,0, ±4^2 


Yes 




/8 0\ 
4 
8 

\0 OJ 




0,0,0,8 


0,8, ±2-^/5 


0,0,±2-\/5 


No 




/ 5 1 -3 1 \ 
1-315 
15 1-3 

V-3 1 5 1 J 




0,0,0,4 


0,4, ±4\/2 


0,0,^:4^ 


Yes 




/ 5 1 -1 3 \ 
3-115 
15 3-1 

V-1 3 5 1 J 




0,0,0,8 


0,8, ±2^ 


0,0,±2V5 


No 



Table 3: Matrices for which neither UECSMTest nor StrongAngleTest are 
applicable. The third and fourth matrices listed are, respectively, unitarily 
equivalent to the first and second matrices. 
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